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Ktedonobacter racemifer corrig. Cavaletti et al. 2007 is the type species of the genus Ktedo- 
nobacter, which in turn is the type genus of the family Ktedonobacteraceae, the type family 
of the order Ktedonobacterales within the class Ktedonobacteria in the phylum 'ChloroflexF. 
Although K. racemifer shares some morphological features with the actinobacteria, it is of 
special interest because it was the first cultivated representative of a deep branching unclassi- 
fied lineage of otherwise uncultivated environmental phylotypes tentatively located within 
the phylum 'Chloroflexi' . The aerobic, filamentous, non-motile, spore-forming Gram-positive 
heterotroph was isolated from soil in Italy. The 13,661,586 bp long non-contiguous finished 
genome consists of ten contigs and is the first reported genome sequence from a member of 
the class Ktedonobacteria. With its 1 1,453 protein-coding and 87 RNA genes, it is the largest 
prokaryotic genome reported so far. It comprises a large number of over-represented COGs, 
particularly genes associated with transposons, causing the genetic redundancy within the 
genome being considerably larger than expected by chance. This work is a part of the Ge- 
nomic Encyclopedia of Bacteria and Archaea project. 



Introduction 

Strain S0SP1-21T [= DSM 44963 = NRRL B-41538) 
is the type strain of the species Ktedonobacter ra- 
cemifer, which is the type species of the monotypic 
genus Ktedonobacter, the type genus of the family 
Ktedonobacteraceae [1]. K racemifer was first de- 
scribed in 2006 [1,2] as an aerobic, non-motile, 
filamentous, mesophilic. Gram-positive hetero- 



troph also capable of growing under microaero- 
philic conditions [1]. The genus name was derived 
from the Greek word ktedon -onos, fiber, and the 
Neo-Latin bacter, a rod, meaning a filamentous 
rod [1]. The species epithet is derived from the 
Latin adjective racemifer, carrying clusters of 
grapes [1]. 
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The original spelling, Ktedobacter racemifer was 
corrected in 2007 on validation according to Rule 
61 and Recommendation 6(7) [2]. Strain SOSPl- 
21^ was originally isolated from a soil sample of a 
black locust wood in Gerenzano, Northern Italy. 
Ten phylogenetically [class level) related strains 
were also isolated from soil samples collected at 
different locations in Northern Italy [1]. Only re- 
cently, a nearest cultivated neighbor, Thermospo- 
rothrix hazakensis, was isolated from hot compost 
in Japan [3]. Here we present a summary classifi- 
cation and a set of features for K. racemifer strain 
S0SP1-21T together with the description of the 
complete genomic sequencing and annotation. 

Classification and features 

Using NCBI BLAST [4], a representative genomic 
16S rRNA sequence of K. racemifer S0SP1-21T was 
compared under default settings [e.g., considering 
only the high-scoring segment pairs [HSPs) from 
the best 250 hits) with the most recent release of 
the Greengenes database [5] and the relative fre- 
quencies of taxa and keywords [reduced to their 
stem [6]) were determined, weighted by BLAST 
scores. The most frequently occurring genus was 
'Ktedobacter' [100.0%) [1 hit in total; this 
represents the original, incorrect spelling of Ktedo- 
nobacter). No hits to sequences with [other) spe- 
cies names were found. [Note that the Greengenes 
database uses the INSDC [= EMBL/NCBI/DDBJ) 
annotation, which is not an authoritative source for 
nomenclature or classification.) The highest- 
scoring environmental sequence was AM180157 
['New Hneage filamentous spore-forming soil iso- 
late SOSP1-30SOSP1-30 str. SOSPl-30'), which 
showed an identity of 99.0% and an HSP coverage 
of 95.2%. The most frequently occurring keywords 
within the labels of environmental samples which 
yielded hits were 'soil' [11.2%), 'prari, tallgrass' 
[4.9%), 'miner, weather' [1.9%), 'new' [1.8%) and 
'filament, lineag, spore-form' [1.6%) [249 hits in 
total). These ke)words reflect some of the ecologi- 
cal properties reported for strain S0SP1-21T in the 
original description [1]. Environmental samples 
which yielded hits of a higher score than the high- 
est scoring species were not found. 

Figure 1 shows the phylogenetic neighborhood of K 
racemifer in a 16S rRNA based tree. The sequences 
of the eight 16S rRNA genes copies in the genome 
differ by up to nine nucleotides from each other and 
by up to five nucleotides from the previously pub- 
lished 16S rRNA sequence [AM180156), which con- 
tains two ambiguous base calls. 



K racemifer strain S0SP1-21T cells are rod- 
shaped, filamentous and grow both vegetative and 
aerial mycelia on solid medium [Figure 2a). The 
large aerial hyphae produce spherical spores that 
cluster together with a grape-like appearance 
[Figure 2b). All other K racemifer strains pro- 
duced rounded spores, although they were ar- 
ranged differently on the aerial hyphae [1]. Fila- 
mentous growth of strain S0SP1-21T also oc- 
curred in submerged cultures, which contained 
the branched mycelia known from actinomycetes 
[1]. S0SP1-21T stains Gram-positive and is not ac- 
id fast [1]. It produces pigments ranging from 
cream to pinkish orange on all media [1]. Although 
essentially aerobic, S0SP1-21T is capable of grow- 
ing under microaerophilic conditions [1]. The op- 
timal growth temperature is 28-33°C [1]. It grows 
well at pH values between 4.8 and 6.8 with an op- 
timum at pH 6 [1]. Salinity up to 10 g per liter does 
not inhibit the growth of the strain [1]. 

Strain S0SP1-21T was capable of hydrolyzing 
starch, casein, gelatin, and [to a lesser extent) kera- 
tin but not cellulose, xylan, or chitin [1]. Strain 
S0SP1-21T was catalase positive and produced H2S 
but could not reduce nitrates [1]. It is sensitive to 5 
ug/ml novobiocin or ramoplanin and to 20 mg/ml 
apramycin and the glycopeptide A40926. 

Chemotaxonomy 

The peptidoglycan of strain S0SP1-21T contains 
ornithine, alanine, glutamic acid, serine, and glycine 
at a molar ratio of approximately 0.7:1.8:1.0:0.8:1.9 
[1]. Serine was identified at the N-terminus of the 
interpeptide bridge [1]. When originally described, 
a detailed peptidoglycan structure had not been 
determined but A-type cross-linkage was suggested 
[1]. The cellular fatty acid pattern of strain SOSPl- 
2 IT was reported to be characterized by an unusual 
high abundance of Cie-.i 2-0H [30%) with other domi- 
nant lipids being branched-chain saturated fatty 
acids /so-Ci7:o [25%), /50-Ci6:o [11.5%) and anteiso- 
Ci7:o [9.6%), as well as Ci6:o 10-Me [7.8%) and Ci6:o 
[6.7%)) [1]. Our own data [DSMZ) did not confirm 
this fatty acid spectrum, but revealed iso-Cie-.o 
[20.1%) as the most frequent fatty acid, followed 
by anteiso-Ci7:o [18.5%), /50-Ci7:o [15.0%), only 
13.1% Ci6:i 2-oH and 11.6% Ci6:o 10-Me. Polar lipids 
consisted of phosphatidylinisitol, phosphatitylgly- 
cerol, diphosphatidylglycerol and an unknown gly- 
cohpid [1]. MK-9[H2) was the only menaquinone 
reported for strain S0SP1-21T [1]. 
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Leptolinea tardivitalis (AB 109438) 
Longilinea arvoryzae (AB243673) 
Levilinea saccharolytica (AB 109439) 



Bellilinea caldifistulae (AB243672) 

Anaerolinea thermolimosa (AB1 09437) 

Anaerolinea thermophila (AB046413) ' 
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Chloroflexus aggregans (CP001337) ** 

0/100 

Chloroflexus aurantiacus (D38365) ** 

Heliothrix oregonensis (L04675) 
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I Herpetosiphon geysericola (AF039293) 

\ 100/100 

' Herpetosiphon aurantiacus (CP000875) ** 



0.03 



Figure 1. Phylogenetic tree highlighting the position of K. racemifer relative to the other type strains within 
the phylum 'Chloroflexi'. The tree was inferred from 1,359 aligned characters [7,8] of the 16S rRNA gene se- 
quence under the maximum likelihood (ML) criterion [9]. Rooting was done initially using the midpoint me- 
thod [10] and then checked for its agreement with the current classification (Table 1). The branches are 
scaled in terms of the expected number of substitutions per site. Numbers above branches are support values 
from 750 ML bootstrap replicates [11] (left) and from 1,000 maximum parsimony bootstrap replicates [12] 
(right) if larger than 60%. Lineages with type strain genome sequencing projects registered in GOLD [13] are 
labeled with one asterisk, those also listed as 'Complete and Published' with two asterisks [14-17] as well as 
CP001337, CP000804, CP000909, CP002084, and AP012029. 
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Table 1. Classification and general features of K. racemifer SOSPl-21^ according to the MIGS recommendations [18] 
and the NamesforLife database [19]. 



MIGS ID Property 



Term 



Evidence code 



MIGS-22 



MIGS-6 

MIGS-15 

MIGS-14 



MIGS-4 

MIGS-5 

MIGS-4.1 

MIGS-4.2 

MIGS-4.3 

MIGS-4.4 



Current classification 



Gram stain 

Cell shape 

Motility 

Sporulation 

Temperature range 

Optimum temperature 

Salinity 

Oxygen requirement 
Carbon source 
Energy metabolism 
Habitat 

Biotic relationship 
Pathogenicity 
Biosafety level 
Isolation 

Geographic location 

Sample collection time 

Latitude 

Longitude 

Depth 

Altitude 



Domain Bacteria J AS [20] 

Phylum Chloroflexi J AS [21 ,22] 

Class Ktedonobacteria J AS [1 -3] 

Order Ktedonobacterales J AS [1,2] 

Family Ktedonobacteraceae J AS [1,2] 

Genus Ktedonobacter J AS [1,2] 

Species Ktedonobacter racemifer J AS [1 ] 

Type strain SOSP1 -21 TAS [1 ] 

positive TAS [1] 

filamentous TAS [1] 

non-motile TAS [1] 

spherical spore-forming TAS [1] 

mesophile TAS [1] 

28-33°C TAS[1] 
NaCi up to 1 0 g/l growth w/o problem, inhibited at 30 g/l TAS [1 ] 

aerobic and microaerophilic TAS [1] 

sugars and peptides TAS [1 ] 

heterotrophic TAS [1] 

soil TAS[1] 

free-living NAS 

none NAS 

1 TAS [23] 

soil from a black locust wood TAS [1 ] 

Gerenzano, Northern Italy TAS [1] 

November 2001 NAS 

45.64 NAS 

9.00 NAS 
not reported 

about 210 m NAS 



Evidence codes - TAS: Traceable Author Statement (i.e., a direct report exists in the literature); NAS: Non-traceable 
Author Statement (i.e., not directly observed for the living, isolated sample, but based on a generally accepted 
property for the species, or anecdotal evidence). These evidence codes are from of the Gene Ontology project [24]. 



Genome sequencing and annotation 
Genome project history 

This organism was selected for sequencing on the 
basis of its phylogenetic position [25], and is part 
of the Genomic Encyclopedia of Bacteria and Arc- 
haea project [26]. The genome project is depo- 
sited in the Genomes OnLine Database [13] and 



the complete genome sequence is deposited in 
GenBank. Sequencing, finishing and annotation 
were performed by the DOE Joint Genome Insti- 
tute QGl). A summary of the project information is 
shown in Table 2. 
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Table 2. Genome sequencing project information 



MIGS ID 


Property 


Term 


MIGS-31 


Finishing quality 


Non-contiguous finished 


MIGS-28 


Libraries used 


Two Sanger 8 kb pMCL200 and fosmid libraries; 
one 454 pyrosequence standard library 


MIGS-29 


Sequencing platforms 


ABI3730, 454 GS FLX 


MIGS-31. 2 


Sequencing coverage 


10.1 X Sanger; 24.6 x pyrosequence 


MIGS-30 


Assemblers 


New^bler version 1 .1 .02.1 5, phrap 


MIGS-32 


Gene calling method 


Prodigal 1 .4, Genemark 4.6b, tRNAScan-SE-1 .23, 
internal 0.81 




INbUC lU 


A / r\r\r\r\r\r\r\r\ 
AUVUUUUUUUUU 




Genbank Date of Release 


June 14, 2010 




GOLD ID 


G1O226I 




NCBI project ID 


27943 




Database: IMG-GEBA 


648276680 


MiGS-13 


Source material identifier 


DSM 44963 




Project relevance 


Tree of Life, GEBA 



Growth conditions and DNA isolation 

K. racemifer S0SP1-21T DSM 44963, was grown in 
DSMZ medium 65 [GYM Streptomyces medium) 
[27] adjusted to pH 6.0, at 28°C. DNA was isolated 
from 0.5-1 g of cell paste using Qiagen Genomic 
500 DNA Kit (Qiagen 10262) following the manu- 
facturer's protocol, with cell lysis protocol 
st/LALMP as described in Wu et al [26]. DNA is 
available through the DNA Bank Network [28]. 

Genome sequencing and assembly 

The genome was sequenced using a combination 
of Sanger and 454 sequencing platforms. All gen- 
eral aspects of library construction and sequenc- 
ing can be found at the JGI website [29]. Pyrose- 
quencing reads were assembled using the Newb- 
ler assembler [Roche). The initial Newbler contigs 
were broken into 14,080 overlapping fragments of 
1,000 bp and entered as pseudo-reads into the 
subsequence assembly. The sequences were as- 
signed quality scores based on Newbler consensus 
q-scores with modifications to account for overlap 
redundancy and to adjust inflated q-scores. A hy- 
brid 454/Sanger assembly was produced using 
parallel phrap (High Performance Software, LLC). 
Possible mis-assemblies were corrected with Dup- 
finisher [30], or transposon bombing of bridging 
clones [Epicentre Biotechnologies, Madison, WI) 
[31]. Some gaps between contigs were closed by 



editing in Consed [32], custom primer walking or 
PGR amplification. A total of 3,354 Sanger finish- 
ing reads and five shatter libraries were produced 
to close gaps, to resolve some repetitive regions, 
and to raise the quality of the finished sequence. 
Illumina reads were also used to correct potential 
base errors and increase consensus quality using a 
software Polisher developed at JGI [33]. The error 
rate of the completed genome sequence is less 
than 1 in 100,000. Together, the combination of 
the Sanger and 454 sequencing platforms pro- 
vided 34.7 X coverage of the genome. The final 
assembly contained 165,050 pyrosequence and 
2,305,667 Illumina reads. 

Genome annotation 

Genes were identified using Prodigal [34] as part 
of the Oak Ridge National Laboratory genome an- 
notation pipeline, followed by a round of manual 
curation using the JGI GenePRIMP pipeline [35]. 
The predicted CDSs were translated and used to 
search the National Center for Biotechnology In- 
formation [NCBI) non-redundant database, Uni- 
Prot, TIGR-Fam, Pfam, PRIAM, KEGG, COG, and In- 
terPro databases. Additional gene prediction anal- 
ysis and functional annotation were performed 
within the Integrated Microbial Genomes - Expert 
Review (IMG-ER) platform [36]. 
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Genome properties 

The non-contiguous finished genome consists of 
ten contigs ranging in size from 1,579 bp to almost 
four Mbp, with five contigs being longer than one 
Mb (1,302,518 bp, 2,713,222 bp, 2,766,182 bp, 
2,916,502 bp, and 3,837,106 bp) and a G+C con- 
tent of 53.8% (Table 3 and Figure 3). Of the 
11,540 genes predicted, 11,453 were protein- 



coding genes, and 87 RNAs; No pseudogenes were 
identified. The majority of the protein-coding 
genes (61.2%) were assigned a putative function 
while the remaining ones were annotated as hypo- 
thetical proteins. The distribution of genes into 
COGs functional categories is presented in Table 4. 



Table 3. Genome Statistics 



Attribute 


Value 


% of Total 


Genome size (bp) 


13,661,586 


100.00% 


DNA coding region (bp) 


10,422,932 


76.29% 


DNA G+C content (bp) 


7,348,426 


53.79% 


Number of contigs 


10 




Extrachromosomal elements 


unknown 




Total genes 


11,540 


100.00% 


RNA genes 


87 


0.75% 


rRNA operons 


8 




Protein-coding genes 


11,453 


99.25% 


Pseudo genes 


0 




Genes with function prediction 


7,065 


61.22% 


Genes in paralog clusters 


4,919 


42.63% 


Genes assigned to COGs 


6,654 


57.66% 


Genes assigned Pfam domains 


7,250 


62.82% 


Genes with signal peptides 


2,660 


23.05% 


Genes with transmembrane helices 


2,581 


22.27% 


CRISPR repeats 


7 
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Table 4. Number of genes associated with the general COG functional categories 



Code 


value 


%aee 


Description 


J 


224 


2.9 


Translation, ribosomal structure and biogenesis 


A 


0 


0.0 


RNA processing and modification 


K 


893 


11.6 


Transcription 


L 


975 


12.6 


Replication, recombination and repair 


B 


3 


0.0 


Chromatin structure and dynamics 


D 


34 


0.4 


Cell cycle control, cell division, chromosome partitioning 


Y 


0 


0.0 


Nuclear structure 


V 


215 


2.8 


Defense mechanisms 


T 


617 


8.0 


Signal transduction mechanisms 


M 


257 


3.3 


Cell wall/membrane/envelope biogenesis 


N 


20 


0.3 


Cell motility 


Z 


0 


0.0 


Cytoskeleton 


W 


0 


0.0 


Extracellular structures 


u 


54 


0.7 


Intracellular trafficking, secretion, and vesicular transport 


o 


195 


2.5 


Posttranslational modification, protein turnover, chaperones 


c 


416 


5.4 


Energy production and conversion 


G 


612 


7.9 


Carbohydrate transport and metabolism 


E 


474 


6.2 


Amino acid transport and metabolism 


tF 


135 


1.8 


Nucleotide transport and metabolism 


H 


264 


3.4 


Coenzyme transport and metabolism 


1 


236 


3.1 


Lipid transport and metabolism 


P 


255 


3.3 


Inorganic ion transport and metabolism 


Q 


217 


2.8 


Secondary metabolites biosynthesis, transport and catabolism 


R 


1,098 


14.4 


General function prediction only 


S 


519 


6.7 


Function unknown 




4,886 


42.3 


Not in COGs 



Insights from the genome sequence 
Genome structure 

With a length of 13;661,586 bp for the ten contigs 
[Table 3) K. racemifer S0SP1-21T has the largest of 
all completely sequenced 1,760 archaeal and bac- 
terial genomes [37] thus far, followed by Soran- 
gium cellulosum, 13.0 Mbp [38], Steptomyces bing- 
chenggensis, 11.9 Mbp [39], Catenulispora acidi- 
phila, 10.5 Mbp [40], and Streptosporangium ro- 
seum, 10.4 Mbp [41]. However, this genome was 
also one of the most difficult to assemble. Figure 4 
shows the unusually high number of identical se- 
quence fragments across the genome, which 
caused the termination of the project as non- 
contiguous finished genome without closure of the 
last ten sequence gaps. 



Comparative genomics 

Lacking an available genome sequence of the clos- 
est relative of K. racemifer, Thermosporothrix ha- 
zakensis [3] [Figure 1), the following comparative 
analyses were done with Sphaerobacter thermo- 
philus [42] and Thermomicrobium roseum [43], the 
closest organisms phylogenetically for which 
there are publically available genome sequences 
[15,16]. 

K. racemifer stands out because of its enormous 
genome size of more than 13 Mbp. The genomes of 
5. thermophilus and T. roseum are significantly 
smaller, 3.9 Mbp and 2.9 Mbp, respectively. Whe- 
reas 5. thermophilus and T. roseum have similar 
G+C-contents of 68% and 64%, respectively, the 
G+C-content of the K. racemifer genome is signifi- 
cantly lower [54%). 
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Figure 4. Screen shot from CROSSMATCH [32] indicating the matches between 
sequences within and across the contigs. CROSSMATCH options were - min- 
match 30 - minscore 60. 



The fraction of shared genes in the three genomes 
is shown in a Venn diagram [Figure 5). The num- 
bers of pairwise shared genes were calculated 
with the phylogenetic profiler function of the IMG- 
ER platform [36]. Homologous genes within the 
genomes were detected with a maximum E-value 
of IQ-^ and a minimum identity of 30%. 

A total of 1,393 genes are shared by the three ge- 
nomeS; referring to the whole genome sizes 39% 
and 48% of the genes in S. thermophilus and 



T. roseum have homologs in the three genomes, in 
the case of K. racemifer only 12% of the genes are 
shared by the other two genomes. The pairwise 
comparison of S. thermophilus and T. roseum re- 
vealed 2,249 genes which are shared by these two 
organisms, referring to the whole genomes 64% of 
the S. thermophilus genes and 79% of the T. ro- 
seum have homologous genes in the respective 
other genome. 
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The genome of K. racemifer encodes an enormous- 
ly high number of transposon-associated genes; its 
annotation revealed 601 genes encoding transpo- 
sases, 151 genes encoding integrases and 107 
genes encoding resolvases. The genes coding 
these enzymes are spread over the whole genome 
with some regions having a higher density than 
others. The extremely high number of transposas- 
es is due to several gene copies that are to a great- 
er or lesser extent similar in their sequences. The 
presence of that many mobile elements may ex- 
plain the unusually high number of identical se- 
quence fragments across the genome and the re- 
sulting difficulties occurring during the genome 
assembly. 

Within the 9,539 unique genes of K. racemifer that 
have no detectable homologs in the genomes of S. 
thermophilus and T. roseum [under the sequence 
similarity thresholds used for the comparison) the 
29 genes encoding xylose isomerases appear to be 
especially noteworthy; for 27 of these isomerase 
genes no homologous genes were detected in the 



other two genomes; only one gene was identified 
in T. roseum, and two in 5. thermophilus. The high 
number of xylose isomerase genes suggests a 
strong utilization of pentoses by K. racemifer. To 
date K. racemifer was not tested regarding xylose 
utilization, but the close relative T. hazakensis is 
able to use xylose as the only carbon source [3]. 
Furthermore, a high number of genes encoding 
proteins responsible for resistance against several 
antibiotics were predicted: 61 bleomycin resis- 
tance proteins and 41 aminoglycoside phospho- 
transferases. 

An estimate of the overall similarity between K. 
racemifer, S. thermophilus and T. roseum, was gen- 
erated with the GGDC Genome-to-Genome Dis- 
tance Calculator [44,45]. This system calculates 
the distances by comparing the genomes to obtain 
HSPs [high-scoring segment pairs) and interfering 
distances from a set of formulas [1, HSP length / 
total length; 2, identities / HSP length; 3, identities 
/ total length). Table 5 shows the results of the 
pairwise comparison between the three genomes. 




Figure 5. Venn diagram depicting the intersections of protein sets (total number of derived protein sequences 
in parentheses) of K. racemifer, S. thermophilus and T. roseum. 
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HSP length / 
total length [%] 


identities / 
HSP length [%] 


identities / 
total length [%] 


K. racemifer 


S. thermophilus 


0.57 


86.4 


0.50 


K. racemifer 


T. rosea m 


0.48 


87.2 


0.42 


T. rosea m 


5. thermophilus 


9.41 


83.1 


7.82 




0 20 40 60 80 100 



Rank 

Figure 6. Relative frequencies of the 100 most frequent COGs in the genome of K. racemifer (blue line) 
compared to their expected frequency as estimated using the broken-stick distribution (red line). Over- 
represented COGs are labeled. 



The pairwise comparison [Table 5) of the ge- 
nomes of K. racemifer with S. thermophilus and T. 
roseum revealed that only 0.57% and 0.48% of the 
average of the genome lengths are covered with 
HSPs. The identity within these HSPs was 86.4% 
and 87.2%, whereas the identity over the whole 
genome was only 0.50% and 0.42%, respectively. 
The comparison of T. roseum with S. thermophilus 
revealed that 9.41% of the average of both ge- 
nome lengths are covered with HSPs, with an 
identity within these HSPs of 83.1%. The identity 
over the whole genome is 7.82%. These results 
show how distant the relationship between K. ra- 
cemifer and S. thermophilus and T. roseum, 



respectively, is, if genome sizes are taken into con- 
sideration. 

In order to quantify the differences in gene redun- 
dancy between the three genomes, as well as to 
determine over-represented genes, we used ap- 
proaches based on entropy and the broken-stick 
distribution, respectively, applied to the set of 
genes from either genome assigned to COGs. 
Shannon's entropy [see, e.g., pp. 214, 243 in [46]) 
H can be used as a measure of disorder for dis- 
crete distributions; it is maximum [Hmax) if all cat- 
egories [COGs in our case) are represented by ex- 
actly one item [gene) and then equal to the loga- 
rithm of the number of items [genes). 
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Thus, one can measure the evenness [non- 
redundancy) within such a distribution as H/Hmax and 
the corresponding redundancy as 1.0 - H/Hmax. The 
broken-stick distribution reflects the relative abun- 
dance of a given number of categories within a ran- 
dom population of items [see, e.g., p. 244 and 410 in 
[46]). Over-represented items [here: COGs) are those 
whose real relative frequencies [here: number of 
genes assigned to this COG relative to the total num- 
ber of genes assigned to COGs) are larger than the 
broken-stick value of the corresponding rank within 
the list of frequencies sorted in decreasing order. 
Moreover, the entropy Hexp of the broken-stick distri- 
bution can be used as an estimate for the expected 
entropy, yielding 1.0 - H/Hexp as an alternative meas- 
ure of redundancy [which becomes negative when 
the evenness is larger than expected by chance). 



The 2,022 genes assigned to 1,300 distinct COGs in 
the genome of T. roseum corresponded to an entro- 
py of 6.912, an expected entropy of 6.748 and, 
hence, a redundancy of 9.20% if measured using 
Hmax and of -2.42% using Hexp, whereas S. thermo- 
philus [2,619 genes assigned to 1,383 COGs) 
yielded an entropy of 6.837 [expected: 6.810) and a 
redundancy of 13.14% with Hmax and -0.39% with 
Hexp- In contrast, the 6,654 genes assigned to 1,731 
distinct COGs in the genome of K. racemifer yielded 
an entropy of only 6.455 [expected: 7.034) and a 
redundancy of 26.67% [using Hmax) and 8.24% [us- 
ing Hexp). That is, in contrast to the other two ge- 
nomes the genes within the genome of K. racemifer 
are distributed less even than expected by chance. 




Rank 

Figure 7. Relative frequencies of the 100 most frequent COGs in the genome of 5. thermophilus 
(blue line) compared to their expected frequency as estimated using the broken-stick distribution 
(red line). Over-represented COGs are labeled. 



Figure 6 compares the relative frequencies of the 
COGs in the genome of K. racemifer compared to 
their expected frequency. More than 80 COGs 
were judged as over-represented by this compari- 
son, considerably more than in the genomes of 5. 
thermophilus [33; Figure 7] and T. roseum [[15]; 
Figure 8). A closer look onto the 20 most over- 



represented COGs in K. racemifer, S. thermophilus 
and T. roseum revealed differences between the 
three organisms. Not surprisingly the genes cod- 
ing transposases [COG0675; by far the most fre- 
quent one), integrases [COG3316) and resolvases 
[COG2452) can be found among the over- 
represented COGs in K. racemifer [Figure 6). 
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Our analyses also showed that genes belonging to 
the category COG3344 are over-represented in the 
genome of K. racemifer. COG3344 represents re- 
tron type reverse transcriptases, which are found 
in group II introns. Group II introns are large cata- 
lytic RNA molecules that act as mobile genetic 
elements [47]. They were first identified in mito- 
chondria and chloroplast genomes, but with the 
increasing number of bacterial genome sequenc- 
ing projects, the number of group II intron se- 
quences in the databases also increased. Dai and 



Zimmerly reported in 2003 that a quarter of the 
sequenced bacterial genomes contain group II in- 
trons [48,49]. By using the IMG-ER platform [36] 
we calculated that approximately one third of the 
2,727 sequenced bacterial genomes contain group 
II introns. In the genome of K. racemifer, 34 genes 
coding reverse transcriptases could be identified, 
all of them having the same domain structure with 
the reverse transcriptase domain followed by a 
maturase-specific domain and the C-terminal 
HNH-endonuclease domain. 




0 20 40 60 80 100 

Rank 

Figure 8. Relative frequencies of the 100 most frequent COGs in the genome of T. roseum (blue line) 
compared to their expected frequency as estimated using the broken-stick distribution (red line). Over- 
represented COGs are labeled. 
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